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The effect of a temporal modulation at three times the critical frequency on a Hopf bifurcation 
is studied in the framework of amplitude equations. We consider a complex Ginzburg-Landau 
equation with an extra quadratic term, resulting from the strong coupling between the external field 
and the unstable modes. We show that, by increasing the intensity of the forcing, one passes from an 
oscillatory regime to an excitable one with three equivalent frequency locked states. In the oscillatory 
regime, topological defects are one-armed phase spirals, while in the excitable regime they correspond 
to three-armed excitable amplitude spirals. Analytical results show that the transition between these 
two regimes occurs at a critical value of the forcing intensity. The transition between phase and 
' amplitude spirals is confirmed by numerical analysis and it might be observed in periodically forced 

' reaction-diffusion systems. 
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PLh ; I- INTRODUCTION 

In many cases, the nucleation of spatio-temporal patterns is associated with continuous symmetry breaking, and 
these patterns are thus very sensitive even to small perturbations or external fields. Perturbations may be induced 
by imperfections of the system itself (e.g. impurities), of the geometrical set-up (e.g. the boundary conditions), of 
the control parameters, etc. In addition, external fields may induce spatial or temporal modulations of the control 
^ or bifurcation parameters. In fact, spatially or temporally modulated systems are very common in nature, and the 
Q^ , effect of external fields on these systems has been studied for a long time. As a way of example, the forcing of a large 
' variety of nonlinear oscillators, from the pendulum to Van der Pol or Duffing oscillators, has led to detailed studies of 
Cp the different temporal behaviors that can be obtained. It has been shown that resonant couplings between the forcing 
and the oscillatory modes may lead to several types of complex dynamical behaviors, including quasi-periodicity, 
frequency lockings, devil's staircases, chaos and intermittency ||f|,||. 
. In equilibrium systems, the importance of spatial modulations has been known for a long time. For example, in 
the case of spatial modulations occurring in equilibrium crystals, such as spin or charge density waves, the constraint 
imposed by the periodic structure of the host lattice leads to the now commonly known commensurate- incommensurate 
phase transitions. The transition from the commensurate phase, where the wavelength of the modulated structure is 
a multiple of the lattice constant, to the incommensurate one occurs via the nucleation of domain walls separating 
^ ' domains which are commensurate with the host lattice 

In nonequilibrium systems, the systematic study of the influence of external fields on pattern forming instabilities 
^ is more recent. It has been first devoted to instabilities leading to spatial patterns. For example, the Lowe-GoUub 

experiment |^ showed that, in the case of the electrohydrodynamic instability of liquid crystals, a spatial modulation 
of the bifurcation parameter may induce dis-commensurations, incommensurate wavelengths and domain walls. The 
similarities with analogous equilibrium phenomena rely on the fact that, close to this instability, the asymptotic 
dynamics is described by the minimization of a potential . 

In the case of self-oscillatory systems, however, original effects occur as a consequence of the nonrelaxational 
character of the dynamics. In particular, for wave bifurcations, unstable standing waves or two-dimensional wave 
patterns may be stabilized by pure spatial or temporal modulations of suitable wavelengths or frequencies [pHll[. 
The case of pure temporal modulations in oscillating extended systems has been considered theoretically | r^ , |I4| and 
in experiments on chemical systems forced by periodic illumination |l5|-[l^. The study by CouUet and Emilsson 
of forced Hopf bifurcations is based on amplitude equations of the scalar Ginzburg-Landau type. They considered 
periodic temporal modulations of frequency uje = {n/m){uJo — i^), where (n/m) is an irreducible integer fraction, loq 
is the critical frequency of the Hopf bifurcation, and v is a small frequency shift. Such forcings break the continuous 
time translation down to discrete time translations, and the corresponding amplitude equations become: 
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dtA - (/i + ii^)A + (1 + ia)W'A - (1 + i(3)A\A\' + 7n^"-\ (1) 

{A stands for the complex conjugate of A). In the non forced case (7,1 — 0) and for zero frequency shift {v = 0) 
this dynamical system is of the relaxational type for a = f3 = 0. In more general situations the system follows a 
nonrelaxational dynamics If the forcing intensity 7„ is sufficiently strong, this dynamics admits asymptotically 
stable uniform steady states, corresponding to frequency locked solutions. There are n different frequency locked 
solutions, which only differ by a phase shift of 27r/n. These solutions are always stable in the large forcing limit for 
n = 1,2,3,4, the so-called strongly resonant cases. In this regime, the dynamics resembles some sort of excitability. 
The locked solutions may undergo various types of instabilities |l3|. One of them is of phase type and occurs 
when 1 + aP is sufficiently negative. In this case, competition between phase instability and forcing leads to the 
formation of stripes or hexagonal patterns, with their associated topological defects. If the forcing is decreased, these 
structures break down through spatio-temporal intermittency p3| . On the other hand, in the phase stable regime, 
frequency locked solutions may undergo a variety of bifurcations when forcing is decreased, leading to oscillation, 
quasi-periodicity or chaos (l9| . 

The equivalence between the different frequency locked states makes possible the formation of stable inhomoge- 
neous structures. These structures are composed of domains of the locked states separated by abrupt interfaces. 
Nonrelaxational dynamics may induce interface motion and, in particular, the formation of n-armed spirals, each arm 
corresponding to a different frequency locked solution. 

These phenomena were studied in great detail by Coullet and Emilsson, for n = 1 and n = 2, in one- and two- 
dimensional systems The case n — 4 has been considered theoretically in jl^. Experimentally, resonant phase 
patterns associated with frequency locking have been described for a periodic forcing of the Bclousov-Zhabotinsky 
chemical reaction for n ~ 2 ||l^ and n = 4 [|8|. For n — 3, the existence of three armed rotating spirals in two- 



dimensional systems is only briefly mentioned in 1 13j . Experimental observations of frequency locking for n = 3 are also 
briefly discussed in ||l^,^. In these experiments patterns with three-phase domains shifted by 27r/3 are observed [ pO[ . 
It is the aim of this paper to study the case of resonant forcing for n = 3 in more detail, and specially the transition 
from phase spirals to amplitude spirals. The interest of this study is threefold. First, it confirms the robustness of the 
Ginzburg-Landau dynamics, which is recovered at low forcing, with all its complexity and its particular sensitivity to 
kinetic coefficients. Second, it presents original dynamical behavior in the excitable regime. This behavior presents 
interesting analogies with Rayleigh-Benard convection in a rotating cell, described by three-mode dynamical models 
pl|-p3[. Third, it could be a useful framework for the interpretation of detailed experiments on the 3:1 resonant 
forcing of a chemical system |^ . 

The pap er is org aniz ed as follows. The dynamical model and its uniform asymptotic solutions are presented in 
section O. Section III is devoted to the description of the dynamics in terms of phase equations. The properties 



and possible development of front and spiral solutions are discussed in section IV. Numerical results, for one- and 
two-dimensional systems, are presented in section ^ and conclusions are drawn in section VI, 



II. UNIFORM SOLUTIONS 

Consider an extended system undergoing a Hopf bifurcation at zero wavenumber, and subjected to a periodic 
temporal modulation of frequency uje = Sujq- Sufficiently close to the bifurcation, its dynamics may be reduced to the 



following complex Ginzburg-Landau equation |24 



dtA = iiA + (1 + ia)\/'^A - (1 + ifi)A\A\^ + 7^^ (2) 

where 7 > (the case 7 < follows by changing A —A) is proportional to the external field intensity. The other 
parameters are standard |l^j2^. We will restrict ourselves to this case of resonant forcing (i/ = 0). A slightly off- 
resonant forcing is known to induce a richer dynamical behavior whose characterization for a spatially extended 
system is beyond the scope of this paper. 

We look now for uniform solutions. By dropping the spatial derivative terms, the corresponding uniform equations 
are, in phase and amplitude variables {A — Ro{t)e^'^°''*^): 

Ro = ^Ro - i?o + 7i?o cos 3$o, 

$0 - -/3i?o-7fiosin3$o. (3) 
If we look at the stationary solutions (fixed points), eqs. give 

il + p^)Rt~i2fi + j^)Rl + ^i^ = 0, (4) 
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from where the amphtudes of the uniform solutions are given by: 



Rl = ^(1 + /3')[2m + 7' ± + 4m7^ - (5) 



Such solutions exist provided that 7 > 7c, with 



7,^-2A*(yTT^-l). (6) 

Note that for /3 = these solutions exist for any nonvanishing forcing. We will consider the case of /? 7^ for which 
these solutions appear at a finite value of the forcing. Once the amplitude is determined by (H), the phase can be 
obtained from the stationary version of (^): 

cos 3$o = , sin 3$o - . (7) 

7-Ro 7 

Each value of Rq — i?_ gives rise to three solutions for the phase $0 which only differ by a phase shift of 2tt/3. 
Hence, for 7 > 7c the system has six uniform solutions: $2 j ^3) = (^1 ; '^'i + 27r/3, $" + 4tt/3) corresponding to 
R- and {^l, $2, ^3) = {^1,^1 + 27r/3, $1 + 47r/3) corresponding to R+. A Hnear stability analysis shows that the 
are always linearly unstable whereas the are stable for \f3\ < \/3. The three $^ solutions are called the frequency 
locked solutions. These become oscillatory unstable {k — 0, iv 0) for \f3\ > VS in the range of forcing amplitudes 
7c ^ 71 < 7 < 72 where 



71 = W7c + TTTTT^ 



4Vl + /32(l - 3/32) _^ 7^^3 _ /32 _^ 3^3^ _ 5 



2(3/32 - 1) 

72 = Vm(1 + /32)/2. (8) 

In the case where the frequency locked solutions are stable, we can show that the system behaves as an excitable 
one: let us construct the nuUclines of the dynamical system (|^), defined as the curves Rq — and $0 — 0, and 
represent them in figure 0, ioi- ■y ~ /3 — I, fi — 0.25, or in figure |2[ for 7 = /3 = 0.01, ^ — 0.25. In both figures, 
it is easy to see that the i?_ states (labeled u) are unstable, while the R+ states (labeled e) are stable for small 
perturbations. However, for perturbations larger than a well defined threshold, the latter are unstable and the system 
makes an excursion in the phase space, before reaching another, equivalent, steady state. It is a form of excitability. 
For 7 < 7c there are no fixed points of (^ and asymptotic solutions correspond to temporal oscillations of the limit 
cycle type. For 7 = the limit cycle is a circle that becomes deformed for < 7 < 7c (see figure ^) . On increasing 
7, the period of the oscillations increases and diverges for 7 ^ 7c. The stability of these oscillatory solutions can be 
better analyzed in the framework of phase equations that we will develop in the next section. 

The transition between locked and oscillatory states is similar to the Andronov-van der Pol bifurcation that 
appears in several types of excitable systems |2q ]. In fact, on decreasing the amplitude of the forcing, stable and 
unstable locked states merge via inverse saddle-node bifurcations which give rise to the birth of limit cycle oscillations. 
Note that we have here three pairs of fixed points merging, instead of just one pair in classical cases. When 7 approaches 
7c from below, the period of the oscillations diverges as (7c — 7)^^^ showing some kind of critical slowing down. On 
the other hand, when 7 > 7c, small perturbations around the stable states decay, while sufficiently large perturbations 
put the system on an heteroclinic trajectory that goes from an stable fixed point to a different, but equivalent fixed 
stable point through the excitability threshold. The intensity of the perturbations leading to an heteroclinic trajectory 
decreases when 7 approaches 7c, and this phenomenon reflects a type of excitable regime. 

All these features can be analytically shown in the limiting case /J, 7 <C ^. In this limit, and taking into account 
that 7c — (3y/JI, the adiabatic elimination of the amplitude in (^ leads to the phase equation 

^0 = -Vm(7c + 7sin3$o), (9) 

from where the excitable stable steady states are given by sin3<I'^ = — 7c/7 and cos3$f > 0, and the three unstable 
steady states satisfy sin3$" = —jc/j and cos3$f < 0, i = 1,2,3. In this case the threshold of the perturbation 
leading to excitability is thus given by: 

A$ = |$5=-$"|~---^. (10) 
' ' * ' 3 3 7 ^ ' 
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Equation (^) describes a relaxational motion of a fictitious test particle in a potential 
where the potential is 

F($o)-$o-:r^cos3$o (12) 
37c 

This potential picture gives a framework to describe the bifurcation at 7 = 7c in the language of phase transitions. 
In the case 7 < 7c the potential is unbounded and has no local minima. As a consequence, the motion is such that 
the phase ^o{t) decreases monotonically and the test particle does not stop in any selected value of the angle. When 
considered modulus 2tt, the trajectory is a periodic one $o(t+r) = $o(t) + 27r with a period T = {2Tr/ ^){-f^-j'^)-'^/'^ 

which diverges as T ~ (— e)^^/^ when j ^ with e = 7/7C ~ 1. As in mechanical problems, it is possible to define 
a probability density function, P(<I>o), for the angle variable as proportional to the time the particle spends in 
the neighborhood of any value of the angle (in other words, the probability is inversely proportional to $0). This 
probability density function P{^o) develops three peaks which sharpen as 7 increases, indicating three preferred 
values for the angle. This shows that the characteristic time around each minimum increases for 7 ^ 7^ and actually 
it diverges as (— e)~^/^. However, for all values of 7 < 7c, the dynamics is such that none of the preferred values for 
the angle is actually selected since the system moves continuously from one to another. In that sense, we can say that 
there is a dynamical restoring of the angular symmetry. 

The situation is quite different for 7 > 7c. In this case, the potential V{^o) develops local minima which, at lower 
order in e are tt/2 + (%/2/3)ei/2, $1 = + 27r/3, = $f + An/S. These minima correspond to the three 

frequency locked solutions. The relaxational motion in the potential is such that for each trajectory, and depending 
on the initial condition, the test particle selects asymptotically one of the local minima and the angular symmetry is 
now broken. The probability density function for that particular trajectory becomes a delta function centered around 
the selected value of the angle. For an ensemble of different initial conditions, the probability density function is a 
sum of three deltas, P($o) = 5 J2^=l ^i^o - $f )■ 

As usual, the transition between the symmetric and the symmetry-broken phases can be characterized by an order 
parameter. For an angular variable, the order parameter should be some periodic function, such as the mean value 
of the sine, (sin(<I>o)) (other periodic functions give similar results). For 7 < 7c it is (sin(<I'o)) = whereas for 7 > 7c 
that average sets into one of three possible values corresponding to the selected phase. Since its value in any phase is 
different from zero at 7 = 7c, the order parameter changes discontinuously at the transition point and, following the 
standard notation, the transition can be classified as a first-order one. In this mean-field description, and according 
to the discussion above, the transition implies a critical slowing down with a characteristic time diverging with a 
exponent 1/2 both below and above the transition point. Equivalently, one can characterize the transition point 
7 = 7c by the fact that the frequency of the motion tends to zero continuously at that point and stays equal to zero 
for 7 > 7c. In an extended system, one could observe coexistence of the three phases at different locations in space, as 
described later in figure ^. Furthermore, in any finite system, there is no true phase transition and the sharp behavior 
predicted by this simple mean-field analysis gets smeared out and the delta-type probability density functions for 
7 > 7c have a finite width. Evidence of this fact is given later from our numerical simulations (see the lower row of 
figure ^. 

III. PHASE APPROXIMATION 

In this section we present several phase equations each one valid in a different region of parameters. As mentioned 
in the previous section, phase equations can be used to analyze the stability of the uniform patterns. 

A. The oscillatory regime 

In the oscillatory regime (7 < 7c), the phase dynamics can be obtained by perturbing the uniform solution 
{Ro{t),^o{t)), and writing R = Ro{t) + p{r,t), $ = <i>o(t -f cf){r,t)). Following Hagan |2^, the adiabatic elimina- 
tion of the amplitude perturbations in the regime /3,7 <C /i leads to the following phase dynamics: 

dt(l>^{l + ap)V^4> + n{V4>f ^ (13) 
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where T is the period of the oscillations, and 



5 JO 2i?,o— 7 cos "Po " JO 2/?,o— 7 cos tPo ' ^' (AA\ 

For 7^0, one recovers the usual Burgers equation 

dt4> = (1 + a/3)V20 + (a - /3)(V^)2 + ... (15) 

with (f) = Pfi(j). 

Hence, in the regime where l + a/3 > 0, stable (phase) spiral waves may be expected, with wavenumber proportional 
to K, and thus depending on the characteristics of the oscillations [ pTfp^ . In this regime, the qualitative behavior and 
interaction between these topological defects should thus be almost insensitive to the forcing p9|-pl|. Furthermore, in 
the regime where 1 + a/? < 0, defect mediated turbulence should also be expected ||3^. In the oscillatory regime, the 
system presents thus qualitatively the same complexity and the same spatio-temporal behaviors than self-oscillating 
systems. Only quantitative aspects are affected by the forcing. 

B. The excitable regime 

In the excitable regime 7 > 7c, and the phase dynamics can be obtained in the limit /3, 7 <C /x, /3 ^ 1 by eliminating 
adiabatically the amplitude of the field. Taking into account that, in this regime, i?^ ~ /i and 7c — we are left 

with the following phase equation: 

dt^ = -^i{ic -f 7 sin 3$) + (1 + a/?)V2$ - (a - /?)(V4')2 + ^^!il±^V''$ (16) 

Besides the homogeneous solutions discussed in the previous section, this equation admits front solutions connecting 
stable states asymptotically at a; = ±cxi. In the case a — P — Q the phase equation is relaxational and the fronts 
connect two states with the same value of the potential and are, therefore, stationary. In the case a = /? 7^ 0, the 
phase equation is still relaxational but now the steady states have different value of the potential and the front moves. 
Moreover, when a ^ (3 there is a purely nonpotential induced front motion. Equation ( p^ ) will be used in the next 
section as the starting point to compute the velocity of the front solution. 

In order to study pattern forming instabilities, we can use ( [l^ ) in the limit of small $. Expanding the sinSfE" up 
to linear order in $, we are led to a damped Kuramoto-Sivashinsky phase equation [p^ . It follows that frequency 
locked solutions are stable for 1 -I- afi > 0. If 1 -I- a(3 < 0, a pattern forming instability of the locked states would 
occur for ]l|[| 

^ > (1 + aPY • ^^'^ 

Since, in this regime, 7 > 7c, a necessary condition for this instability is thus 



(l-Ha/3)^>72(v/^+^-l)«*(l-t-/3')^ (18) 

and this condition cannot be realized in the {1 -t- a[3 < 0} domain. Therefore, the frequency locked solutions are 
stable so that pattern forming instabilities are ruled out within the phase approximation. It is also possible to prove 
that the locked solutions are always stable in the limit of large forcings . 



^This corrects the missprint of reference llSl in e and fco after eq. (28). 
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IV. FRONTS AND SPIRALS 



For 7 > 7c, the forced Ginzburg-Landau equation possesses three equivalent excitable steady states. The excitability 
mechanism described in the previous section provides a natural way of building fronts between these steady states. 
Despite the equivalence of the fixed points, such fronts are expected to move, as a result of the nonpotential character 
of the dynamics. 



A. One-dimensional systems 



Consider a front solution of eq. (16), e.g. <i>i2(a; — vt), joining the states 1 (at a; ^ — oo) and 2 (at x +oo), such 
that $2 > ^1- Its velocity may be computed along the standard procedures, and is such that (at leading order in 
perturbation) 



V 



27r7cV7^ + 3(a-/?)/+^(a,<i>i2)3 



(19) 



12 



Hence, for a > f3, the fronts $12, <&23 and $31 move to the right, while the fronts $21, <&i3 and $32 move to the left. 
Hence, any domain of one steady state, embedded in a domain of another one, either expands or shrinks, leaving the 
system in one steady state (domains of 2 embedded into 1, 3 into 2 and 1 into 3 shrink while domains of 1 into 2, 2 
into 3 and 3 into 1 expand). However, a succession (from left to right) of domains with states in the order 1, 2, 3, 1, 
etc. moves as a whole to the right. When it is in the order 1, 3, 2, 1, etc., it moves as a whole to the left (see figure 



B. Two-dimensional systems 



In two-dimensional systems, straight linear fronts have the same behavior as in one-dimensional systems. Further- 
more, sets of two inclined fronts separating domains with different steady states also move away or annihilate, leaving 
the system in one steady state only (see figure |^) . 

New phenomena may arise when the three steady states coexist in the system. In this case, three fronts, which 
separate the respective domains, coalesce in one point (a vertex). The three fronts are expected to rotate around this 
point. The result is a rotating spiral whose angular velocity increases with the forcing amplitude. Spirals corresponding 
to sequences of states in the order 1^2— >3orl— >3— >2 around the center have opposite senses of rotation. 
Isolated vertices remain immobile, but non-isolated ones have a dynamical evolution induced by mutual interactions, 
which may even lead to the annihilation of counter-rotating spirals. This dynamical behavior is illustrated by the 
results of the numerical analysis presented in the next section. 

The situation is similar to that observed in system with competing fields. For example, in the context of fluid 
dynamics, a three mode model has been proposed to study Rayleigh-Benard convection in a rotating cell |2^j2^. The 
fields represent the amplitudes of three set of convection rolls oriented 60° at each other. In the two-dimensional 
system, vertices may form when the three different types of roll domains meet at one point (notice that this is not 
possible in one dimensional systems [Q). Then, the nonpotential dynamics induces the rotation of the interfaces 
around the vertices preventing the system from coarsening. At long time scales, the vertices diffuse throughout the 
system. 



V. NUMERICAL RESULTS 



In this section, we present numerical results in two spatial dimensions which illustrate the various dynamical regimes 
described in the preceding sections. We have solved numerically the forced CGLE in two spatial dimensions by using 
a pseudospectral method with periodic boundary conditions. We discretize the system in a square mesh of 256 x 256 
points. Cases within and beyond the validity of the phase approximation and with 1 -I- a/3 > and 1 -I- a/3 < are 
considered. In all cases the parameter /i is taken fixed fj, = 1. In table | we summarize the parameters chosen for each 
of the cases studied. 

We start our discussion considering cases in which the phase approximation is valid (/?, 7 ^ fi). We first consider 
a case below the Benjamin-Fair (BF) line (1 + a/3 > 0) choosing parameter values {a — 2, P — —0.2; = 0.2}. 
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This case corresponds to the frozen states regime in the phase diagram of the CGLE with no forcing |3^. In figure 
^ we show the modulus and the phase of the complex field A for values of 7 corresponding to no forcing, oscillatory, 
7 ~ 7c and excitable regimes. As expected [Q, spiral defects surrounded by shock lines occur when there is no 
forcing. When the strength of the forcing 7 is increased, but still being below the critical value jc (oscillatory regime), 
the phase dynamics does not change significantly. However, amplitude spirals appear in the modulus of the field. 
The splitting of the phase into three locked states is observed approximately at the predicted theoretical value of 
the forcing 7c. For a value of the forcing parameter slightly greater than 7c, we observe annihilation of vertices 
until an homogeneous state is reached. For values of 7 close to 7c the motion of the walls can be very slow with 
patterns that might look stationary in short time scales of observation. For larger values of the forcing parameter, 
the nonrelaxational dynamics is able to stop vertex annihilation and therefore the coarsening process. The system 
remains in a self-sustained dynamical state dominated by three-armed rotating spirals. 

Motivated by the way in which experimental data for 2:1 or 4:1 resonant forcing is presented in p7|jlq] , we show 
in figure ^ a phase portrait in the complex plane of A, obtained from the snapshots shown in figure^. We also 
include the corresponding histograms for the values of the phase of A. This representation of data clearly displays 
the transition from the oscillatory to the excitable regime. The phase portrait gives a demonstration, for the spatially 
extended system, of the deformation of the limit cycle described in figure |^. Furthermore, it enlightens the transition 
from one-armed spirals, which are generic defects of unforced oscillations, to three armed spirals, which are generic 
defects of the locked states. Effectively, at 7 = 0, one observes typical phase spirals associated with the monotonic 
phase variations. In this case, there are no fronts separating different states of the system. When 7 increases, the 
phase variable spends an increasing amount of time around the precursors of the locked states, as reflected in the 
histograms of Fig. |^. This manifests itself in the appearance of fuzzy three armed patterns around vertices of the 
modulus of the field A (see for example Fig.^ . For 7 slightly larger than 7c the system has undergone the transition 
form oscillatory to locked behavior and one may observe spirals generated by fronts between different stable steady 
states of the system. This transition is in fact the spatial unfolding of the evolution of the phase portrait displayed in 
Fig.0 for increasing 7, which in turn corresponds to the mean field description of the transition given after eq. (|2|). 

Above the BF line {l + a(3 < 0), we choose a = 5.5 and keep the rest of parameters as before. The most noteworthy 
difference with the previous case is the existence of asyniptotic frozen states for all values of 7 within the oscillatory 
regime, even close to the critical forcing 7c (see figure ||). Below 7c, we observe frozen targets while close to the 
transition (7 « 7c) the frozen patterns hold three locked phase states but without vertices. The particular geometry 
of these patterns is due to the emergence of locked phase states. As expected, large enough values of the forcing 
parameter give rise to a time-dependent dynamics with three-armed spirals rotating around vertices. 

Beyond the validity of the phase approximation, different phenomena may occur. In particular, pattern forming 
instabilities may take place for small and moderate values of the forcing parameter above its critical value 7c. In 
the case above the BF line with parameters {a = 2, p = —0.76; 7c = 0.72} (phase turbulence regime in the absence 
of external forcing) , oscillating targets that coexist with vertices are observed close to the transition (see figure |^) . 
Pulses in the modulus of A form at the center of the targets concentric rings with an amplitude which decays in space 
as they propagate away from the center, while the phase oscillates periodically between — tt and tt. 

We also considered, beyond the validity of the phase approximation, a case below the BF line. Our choice of 
parameters {a = 0, /3 = —1.8; 7c = 1.45} corresponds to a defect turbulence regime at 7 = 0. Since 1 -I- a(3 > 0, well- 
developed spirals can be observed (see figure [To| ) . The modulus of the field is characterized for 7 < 7c by amplitude 
spirals that rotate around defects, whereas the phase shows a behaviour similar to the one for 7 = 0. These three 
armed amplitude spirals become well developed for 7 > 7c corresponding to the emergence of the three possible 
homogeneous phase values. On the other hand, since |/3| > a/S, and according to the discussion of section y, there 
exists a range of values of the forcing par ameter for which the locked solutions present a oscillatory instability at 
zero wavenumber. This is seen in figure This instability is observed after the annihilation of two counter-rotating 
defects. In the squared region identified in the figure we observe the development of an oscillating target corresponding 
to the homogeneous oscillation. ^From figure [Ll| d onwards the oscillating regions shrinks and disappears under the 
invasion of neighboring spirals. 

It is important to emphasize that for 7 3> 7c the asymptotic state is essentially the same regardless of the different 
dynamical regimes that exist for 7 = and different parameter values. The phase is locked to either of the three values 
predicted theoretically, with interfaces between these three locked states rotating around vertices. This rotation, which 
is due to the underlying nonrelaxational dynamics, inhibits the coarsening process which would take place through 
vertex annihilation. The vertices are essentially pinned and the resulting pattern is, on the average, time periodic at 
relative short time scales. When the phase approximation is valid, the locked phase states are seen to be stable, but 
excitable spirals may be absent near the transition for system parameters such that 1 -I- a/3 < 0. Beyond the validity 
of the description based on the phase approximation, instabilities of the homogeneous phase states may take place 
giving rise to complex patterns. 
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VI. SUMMARY 



Temporal forcing of nonlinear self-oscillating extended systems strongly couple with with the unstable modes asso- 
ciated with a Hopf instability. Such forcings modify the character of the bifurcation and the resulting spatio-temporal 
patterns, as observed in different resonant regimes of a forced reaction-diffusion system ]l5|-p^. In this paper we have 
studied the particular case of a resonant temporal modulation at three times (n — 3) the critical frequency of the 
Hopf bifurcation. 

For forcing amplitudes 7 below a critical value 7c, the system is in an oscillatory regime, where the spatio-temporal 
behavior strongly depends on the parameters of the associated Ginzburg-Landau equation. Uniform solutions corre- 
spond to temporal oscillations of the limit cycle type and topological defects correspond to one-armed phase spirals. 
When 7 approaches 7c from below, the period of the limit cycle diverges as (7c — 7)^^^ showing some type of criti- 
cal slowing down. For forcing amplitudes above the critical one, the system is in a phase locked regime with three 
equivalent steady states. The bifurcation occurring at 7 = 7c can be described, in a mena field approximation, as a 
first order phase transition between a motion in which the temporal average of the phase is zero (7 < 7c) to one with 
a nonzero average phase (7 > 7c). We have given a phase approximation description of both regimes for a spatially 
extended system. 

Like in the n — \ and n = 2 cases of strongly resonant forcings, a form of excitability may also be observed. 
However, contrary to the n — 1 and n — 2 cases, no pattern forming instability of the frequency locked states occurs, 
in this case, for parameter values for which the phase approximation is valid. Due to the nonrelaxational character 
of the dynamics, fronts between equivalent steady states move. The result is that, when the three equivalent steady 
states coexist in the system, three armed rotating spirals are generated around vertices where the fronts separating 
each domain meet. Hence, we predict a transition from one-armed phase spirals to three-armed excitable amplitude 
spirals, which occurs when the forcing amplitude passes through a critical value 7c. We have confirmed and described 
this transition by numerical analysis of the corresponding Complex Ginzburg Landau equation for different parameter 
values which correspond to qualitatively different regimes of the phase diagram of this equation when there is no 
forcing. 

There are a number of analogies of our results for 7 > 7c and the spatiotemporal patterns observed in a model of 
rotating Rayleigh-Benard convection ||2^ . In the latter case the domains correspond to sets of paralell convection rolls 
with a certain orientation and the vertices to points at which fronts separating domains of three preferred orientatons 
meet. As in the case studied here, the rotation of interfaces around vertices, due to nonrelaxational dynamics, produce 
rotating three armed spirals. 
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excitability threshold 




1 2 3 4 5 6 

FIG. 1. NuUclines and fixed points of the forced Ginzburg-Landau equation for 7 = /3 = 1, /i = 0.25 in the R, (j) plane. Black 

arrows indicate the dynamical flow along the nuUcline R = 0, while grey arrows indicate the dynamical flow along the nuUcline 
$ = 




excitability threshold 




excitability 



FIG. 2. NuUcline i? = in the complex A plane. Along this nullcline the fixed points u and e are determined by $ = 0. 
Parameter values are "f = l3 = 0.01, ^ = 0.25. Note that for these parameter values R+ ~ R-. Arrows indicate the dynamical 
flow, with long arrows indicating the excitable excursions associated with an heteroclinic trajectory. 
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FIG. 3. Uniform solutions of eq. (g) for several values of the forcing parameter 7. System parameters are /i = 1, /3 = —2.0 
(so 7c = 1.57). Note that for 7 < 7c uniform solutions are of the limit cycle type whereas for 7 > 7c they become fixed points. 
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FIG. 4. Plot of the phase field in the excitable regime in Id. Note the existence of three homogeneous phase states. The 
arrows indicate the direction of motion of the several fronts. Parameter values are = 1, a = 2, /? = —0.2, 7 = 0.5. 
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FIG. 5. Examples of motion of pairs of inclined fronts separating domains with equivalent steady states of the forced 
Ginzburg-Landau equation for 7 > 7c. 





a 


/? 


7c 


Phase approx. 
valid? 


l + a(3 


Regime (7 0) 


Figure 


1.0 
1.0 


2.0 
5.5 


-0.20 
-0.20 


0.20 
0.20 


Yes 
Yes 


> 
< 


Frozen states 
Frozen states 


' — ^ 


1.0 
1.0 


2.0 
0.0 


-0.76 
-1.80 


0.72 
1.45 


No 
No 


< 
> 


Phase turbulence 
Defect turbulence 


' — ^ 



TABLE I. Parameters of the various cases discussed in section 
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FIG. 6. Modulus and phase of the complex field A in the cases 7 = (no forcing), 7 < 7c (oscillatory), 7 ~ 7c and 7 > 7c 
(excitable). Parameter values are /i = 1, a = 2, /3 = —0.2 (so that 7c ^ 0.2), and 7 = 0.1 (0.25) for the oscillatory (excitable) 
case. 
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FIG. 8. Same as in figure |6[ Parameter values are = 1, a = 5.5, /3 — —0.2 (7c — 0.2), and 7 = 0.1 (0.25) for the oscillatory 
(excitable) case. 
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FIG. 9. Same as in figure |6[ Parameter values are /i = 1, a — 2, (3 — —0.76 (7c — 0.72), and 7 — 0.5 (1.5) for tlie oscillatory 
(excitable) case. 
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FIG. 10. Same as in figure ^. Parameter values are = 1, a = 0, /3 = —1.8 (7c — 1.45), and 7 = 1 (1.6) for tlie oscillatory 
(excitable) case. 
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